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Abstract 

We investigate the linear and nonlinear electric polarizabilities of small 
AI4M4 (M=Li, Na and K) clusters. Quantum chemical calculations reveal 
that these compounds exhibit an exceptionally high magnitude of linear and 
nonlinear optical (NLO) coefficients which are orders of magnitude higher 
than the conventional 7r-conjugated systems of similar sizes. We attribute 
such phenomenal increase to non-centrosymmetricity incorporated in the sys- 
tems by the alkali atoms surrounding the ring leading to charge transfer with 
small optical gap and low bond length alternation (BLA). Such a low mag- 
nitude of the BLA from a different origin, suggests the possibility that these 
clusters are aromatic in character and along with the large NLO coefficients 
they appear to be better candidates for next generation NLO fabrication 
devices. 
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The development of materials with large nonlinear optical (NLO) prop- 
erties is a key to controlling the propagation of light by optical means. In 
particular, the response of the materials to the application of the electric field 
has found tremendous applications in designing materials for NLO devicesfp. 
These devices are being used in numerous applications, from lasers to optical 
switches and optoelectronics. The NLO properties of organic 7r-conjugated 
materials have been studied in great details in the last few decades [3 El • 
The second and third order non-linear optical properties, /3 and 7 for the 
TT-conjugated polymers increase with the conjugation length (L) roughly as 

and respectively j3]. Therefore, the general strategy to model NLO 
materials has been to increase the conjugation length. However, there ex- 
ist an upper limit for every off-resonant susceptibilities ji^. Alternatives to 
these TT-conjugated compounds are yet to be explored theoretically in a de- 
tailed fashion. But, with the gaining popularity of various ab-initio level 
methods[n], there has been a tremendous impetus in investigating the struc- 
ture and electronic properties of both homogeneous and heterogeneous small 
clusters in recent years j3IH]. 

Small AI4 rings like AI4M4 and their anions AI4M3 , M=alkali metals, 
have been a subject of current interest [HIIIIII because of their unique charac- 
teristics and close structural resemblance with the C4H4. However, although 
C4H4 is an anti-aromatic species, these Al4-clusters are recently reported to 
be a aromatic Thus, it would be interesting to ask whether these rings 
are better polarized than their organic counterpart; whether the structural 
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characteristics has any role in their polarization response functions. Organic 
TT-conjugated systems are stabilized due to vr-electron delocalizations, while 
the inorganic metal complexes reduce their energy through strong charge 
transfer. There have been no previous efforts to study in details the NLO 
properties of these all-metal clusters. We describe in the following that these 
metal clusters offer a unique polarization response due to their ionic charac- 
ter, contrary to conventional 7r-conjugated systems, leading to large optical 
coefficients. 

We begin our calculations by optimizing the ground state geometries of 
the Al4-clusters (Al4Li4, Al4Na4 and AI4K4). All the optimizations have 
been done using the 6-31G(d,p) basis set. Electron correlation has been 
included according to the DFT method using Becke's three parameter hybrid 
formalism and the Lee- Yang-Parr functionals (B3LYP) J2] available in the 
GAMESSpSj electronic structure set of codes. Since, we want to compare 
the optical properties of these small four-membered rings with their organic 
analogue C4H4, we start with a planar initial geometry for the optimizations. 
We have varied the level of basis set from 6-31G(d,p) to 6-311G+(d) to 
ensure that these geometries correspond to the minima in the potential energy 
surface. The final geometries indeed remain independent of the selection of 
the basis set. Contrary to that of C4H4 having a rectangular ring, these AI4- 
clusters are found to have a rhomboidal structure with four Al atoms forming 
a rhombus and the four alkali atoms around the four Al-Al bonds forming 
four Al-M (Li, Na, K)-A1 triangles. One of the diagonals of the AI4 ring are 
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also connected by Al-Al bond. The equilibrium geometries are shown in Fig. 
1 [1(a), 1(b) and 1(c)]. 

While Al4Li4 and AI4K4 have a planar structure {D2h), Al4Na4 has a 
distorted structure, with the four Na atoms arranged in a nonplanar geometry 
around the planar ring (the Na atoms are distorted by 13 degrees from the 
plane of AI4 ring). This can be understood by considering the increase in 
size of the alkah ions and the distances of the ions from the Al4-ring. With 
the progressive increase in the ionic radii of counterion, Li to K(Li=0.68 A, 
Na=0.97 A and K = 1.33 A), the structures are expected to be distorted 
and the four alkali atoms should arrange in a non-centrosymmetric geometry 
around the AI4 ring to minimize steric repulsion. But, the average Al-M 
distance increases while going from Al4Li4 (2.65 A) to Al4Na4 (3.00 A) to 
AI4K4 (3.35 A). Although the ionic radii of K ion is more than that of Li and 
Na, in AI4K4, the four K ions are far separated from the AI4 ring, allowing 
a planar structure. For Al4Na4, both the ionic radius of Na and the average 
Al-M distance fall in between Al4Li4 and AI4K4 and thereby minimizes the 
steric repulsion through distortion. 

Also, very close in energy to these planar rhomboidal structures for these 
Al4-clusters are the capped octahedron structures for the Al4Li4 [2(a)] and 
Al4Na4 [2(b)] (with C2h symmetry) and a distorted tricycle-like structure for 
AI4K4 [2(c)]. These geometries are shown in Fig.2. At the footnote of each 
structure, the corresponding ground state energies are given. It has however 
not escaped our attention that previous works on alkali derivatives of AI4- 
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clusters have predicted more than one unique structures for these systems jE]. 
This calls for a study to elucidate whether the optical properties for these AI4- 
clusters for different geometries are substantially different or are very similar. 
Hence, both the geometries for each cluster were considered for computing 
the optical response functions. 

These geometries were used to compute the SCF MO energies and then 
the spectroscopic properties using the Zerner's INDO method^^]. We have 
varied the levels of CI calculations, with singles(SCI) and multi-reference dou- 
bles CI (MRDCl), to obtain a reliable estimate of the second order optical 
response. The later method is particularly important since it includes cor- 
relation effects substantially. The MRDCI approach adopted here has been 
extensively used in earlier works, and was found to provide excitation ener- 
gies and dipole matrix elements in good agreement with experiment |16t IT7j . 
As reference determinants, we have chosen those determinants which are 
dominant in the description of the ground state and the lowest one-photon 
excited states^^. We report the MRDCI results with 4 reference deter- 
minants including the Hartree-Fock HP ground state. For each reference 
determinant, we use 5 occupied and 5 unoccupied molecular orbitals to con- 
struct a CI space with configuration dimension of 800 to 900. To calculate 
NLO properties, we use correction vector method, which implicitly assumes 
all the excitations to be approximated by a correction vector J^l- Given the 
Hamiltonian matrix, the ground state wave function and the dipole matrix, 
all in CI basis, it is straightforward to compute the dynamic nonlinear optic 
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coefficients using either the first order or the second order correction vectors. 
Details of this method have been published in a number of papersj^Hl 1^122] • 

Table 1. shows the bond-length alternation (BLA), Ar, the optical gap 
and the average Mulliken charge on the AI4 ring for all the geometries. The 
Ar is defined as the average difference between the bond lengths of two 
consecutive bonds in the Al^-img and the optical gap is calculated as the 
energy difference between the geometry relaxed ground state and the lowest 
optically allowed state with substantial oscillator strength. This corresponds 
to the vertical absorption gap. To directly compare the efficiency of these 
AU-clusters with the conventional n conjugated systems, we calculate the 
optical properties of the 1,3-cyclobutadiene (C4H4) and benzene (CgHg) at 
the same level of theory. For the Al4-clusters, there is a substantial amount 
of charge transfer from the alkali atoms to the Al atoms (negative charge), 
making them act as donor and acceptor respectively. Such a charge transfer 
induces polarization in the ground state structure and reduces the optical 
gap. On the other hand, the C-H bond being perfectly covalent, there is 
almost no charge transfer in case of C4H4 and CeHg and thus have a large 
optical gap due to finite size molecular architecture. 

Charge transfer stabilizes the system with very small changes in the 
bond lengths. The chemical hardness, i], defined as, one half of (ionization 
potential-electron affinity), decreases as one moves from Li to K. More specif- 
ically, the r] for Li, Na and K are 2.39eV, 2.30eV and 1.92eV respectively|23j. 
So, the extent of charge transfer from the alkali atom to the Al4-ring should 
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increase with the decrease in the chemical hardness of the alkah atoms which 
is evident from Table 1. From Al4Li4 to Al4Na4 to AI4K4, the Mulliken 
charge on the Al4-ring increases leading to decrease in the BLA along the 
series with the exception of Al4Na4 [1(b) and 2(b)] which has much lower 
Mulliken charge on the Al4-ring. For, the Al4Na4 [1(b)] as mentioned above, 
there is a substantial distortion of the Na atoms from the Al4-ring. For the 
Al4Na4, with C2/1 symmetry [2(b)] even though there is no distortion, the 
large Al-Na distance reduces the ionicity of the bond. As a result the extent 
of charge transfer is lesser for Al4Na4. 

Benzene is aromatic with Ar = 0. Although BLA can not be regarded 
as the sole measuring index of aromatic character, the small BLA found for 
the Al4-clusters (together with large BLA for the anti-aromatic C4H4) tend 
to suggest that the Al4-clusters are more like aromatic but most certainly 
not anti-aromatic species, as has been proposed recently ^HI- The distorted 
AI4K4 structure [2(c)] is very interesting. The A^-ring is distorted from 
planarity by 9.5 degrees. Such a distortion arises to minimize the steric 
repulsion in accommodating four bulky K atoms on a plane, very similar 
to that of cyclooctatetraene, which undergoes a distortion from planar to 
tub-shaped geometry to minimize the ring-strain j21]. Thus, this structure 
of AI4K4 is neither aromatic nor anti-aromatic but can be considered to be 
non- aromatic just like cyclooctatetraene. This is supported by the energies 
for the structures for AI4K4 [1(c) and 2(c)]. The distorted structure is more 
stable than the planar structure. Thus, the steric repulsion for the four large 
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K atoms overwhelms the stabihty of the planar undistorted aromatic AI4K4 
making the AI4K4 cluster non- aromatic. 

In Table 2, the magnitudes of the ground state dipole moment, fic, the 
linear (a), and nonlinear (/3 and 7) polarizabilities for the clusters are reported 
from the ZINDO calculations. Note that we report the magnitudes for the 
tumbling averaged d, j3 and 7, defined asj2n] 

^ = ^ 

7 = T7ll(27iw + 7yii) (1) 

where the sums are over the coordinates x,y,z = x,y,z) and P* refers 
to the conjugate of /3j vector. 

The ground state dipole moment fie and P are zero for the C4H4 and 
CgHg due to its perfect centrosymmetric geometry, although, the a and 7 
have finite values. For the A^-clusters with the progressive increase in the 
ionic radii of counterion, the ground state dipole moment increases. Thus, 
while Al4Li4 has no ground state dipole moment, Al4Na4 and AI4K4 have 
substantial ground state dipole moment (particularly 1(b) and 2(c) due to 
their non-centrosymmetric structures discussed above). For AI4K4, while the 
rhomboidal geometry has a very low ground state dipole moment but the 
distorted tricycle-like structure has a very high dipole moment. Thus, due to 
the out-of-plane charge transfer, the dipole matrix elements are also larger, 
resulting in particularly large value for jS. For Al4Li4 and Al4Na4 [2(b)] the 
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polarization is in the excited state as the ground state dipole moment is zero. 
However, it is not the case for the insulating C4H4 and CeHg which have zero 
polarization both in the ground and the optical excited states. Thus, (3 is 
zero for C4H4 and CgHg. 

The optically active states are the low-energy states of these metallic 
clusters and the lowest optical gap is about 0.07 au for Al4-clusters compared 
to 0.25 au for the C4H4 and CgHg. Since the optical coefficients are inversely 
proportional to the optical gaps and proportional to the dipolar matrices, a 
large optical gap implies low magnitudes for the optical coefficients. C4H4 has 
the highest magnitude of BLA and optical gap and the least charge transfer 
on the ring structure, thereby smallest magnitude of 7. On the otherhand, 
although BLA is zero for CeHg due to complete vr-electron delocalization, 
there is no charge transfer in the finite molecular structure leading to large 
optical gap and weak polarization. Consequently, 7 is very less also for CgHg. 

In contrast, the optical coefficients in general are quite large for the AI4- 
clusters. For example, the 7 for the Al4-clusters are roughly 10^ times more 
than that for C4H4 and CeHg. This is because the 7 is a third order prop- 
erty with 4 dipolar matrices in the numerator and 3 optical gaps in the 
denominator |26j . The 7 for the Al4-clusters increases with the increase in 
the polarization of the Al-M bonds and follows the trend: 7 of Al4Li4 < 7 
of Al4Na4 < 7 of AI4K4 (same trend as r]). But, the distorted structures for 
Al4Na4 and AI4K4 [1(b) and 2(c)] have less 7 due to less polarization of the 
Al-M bonds. For C4H4 (Ar= 0.245A), CeHg (Ar= O.OOA) and linear chain. 
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(-CH=CH-)„, n=3 (Ar= O.lA) the 7 are 2.21, 2.63 and 192.85 (all in au) 
per CH bond, respectively. This is in agreement with previous findings that 
the magnitude of 7 varies nonlinearly with Ar and its maxima occurs at an 
optimal Ar 7^ [27j . However, for the Al4-clusters, because of strong charge- 
transfer, the 7 are less sensitive to variation in Ar and even for Ar=O.OA 
(perfect square AI4 ring) and Ar=0.245A (rectangular AI4 ring like C4H4) 
the 7 are similarly high, as found for the optimized structures. Such charge- 
transfer induced large NLO is evident for even bulk materials like CsLiBeOio 
(CLBO) showing significant 4*'^ and 5*^ harmonic generation 

The maximum possible value for the off-resonance 7(Al4M4) calculated 
using the Kuzyk's simple two-state model ^ is about 2825 times more than 
that for 7(04134). This is in very good agreement with our MRDCI calcu- 
lations at a low frequency (0.001 au), which predict 7(Al4M4)/7(C4H4) ^ 

To compare and contrast these clusters with their organic counterparts, 
we calculate the NLO properties of the well-known 7r-conjugated systems, the 
trans-polyacetylene chain, (-CH=CH-)„, by varying the number of spacers, 
n, from n=l to 6, and thereby extending the length of conjugation from 2.65 
au to 29.11 au. The geometries were optimized by the same method as men- 
tioned above. The linear and nonlinear polarizations are calculated at the 
same frequency (0.001 au). Our calculated values for the optical properties 
compare fairly well in trends with the experimental results that the linear 
(a) and nonlinear (7) optical properties increase steadily with the increase in 
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the conjugation- length of the chain (see Table 2 and Table 3). For example, 
for ethylene, ^expt= 1504.9 au, for butadiene (n=2), ^expt= 4566.4 au and for 
hexatriene (n=3), ^expt= 14950.1 au I2H]. Note that our calculations are done 
at a lower frequency compared to the laser frequency used in the experiment 
(0.066 au). However, the magnitudes of all the polarization quantities are 
much higher for the charge-transfer complex (AI4M4 clusters) compared to 
the conventional vr conjugated chains with comparable conjugation length. 
Only when there are very large number of spacers (n=5-6), that the mag- 
nitudes become comparable to the much smaller Al4-clusters. For example, 
7[(-CH=CH-)6] ^ gamma(Al4M4). {-CH = CH-)q has 12 atoms in con- 
jugation while AI4M4 has only 4. So, as a thumb-rule, one can state that 
gamma for the Al-atoms in the charge-transfer ring scale three times that for 
pi-conjugated organic materials. 

As discussed, the large NLO properties for AI4M4 is due to the charge 
transfer from the alkali metals to the AI4 ring. It will be thus of interest to 
compare these hetero atomic all-metal clusters with alkylated organic com- 
pounds such as lithiated benzene or organolithium and organosodium deriva- 
tives like C8H6Li2 and C8H6Na2. These alkylated organic compounds also 
exhibit larger NLO coefficients. For example, 7(C6Li6)/7(C6H6)=7.3 x 10^ 
jSm. Similarly, 7(C8H6Li2)/7(C8H6)=5.5 and 7(C8H6Na2)/7(C8H6)=20 
But, NLO responses for the AI4M4 is much higher than these alkylated or- 
ganic compounds compared to pure organic materials. Also, there have been 
previous efforts to calculate the NLO coefficients in inorganic clusters like 
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GaN, GaP and GaAs [22]. These systems have higher gap than the AI4M4 
clusters so the NLO coefficients are smaller. 

To conclude, our theoretical study shows that the small four membered 
Al4-clusters functionalized with various metal cations provide an innovative 
route for selection of materials with very high nonlinear optical properties. 
We believe that our study will motivate experiments on these Al-clusters. 
Some of these compounds have already been well characterized from stable 
allovsfTU] but laser evaporation is not sufficient to stabilize these materials for 
NLO experiment. One idea will be to stabilize these clusters in a sandwitch 
type geometry by incorporating a suitable transition metal ionjSSj or solvents. 
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Table 1: The bond length alternation, Ar (in ^4), Optical Gap (in au) and 
the average MuUiken charge {Aq) on the ring for the clusters from ZINDO 
calculations. 



Molecule 


Ar 


Gap 


Aq 


Al4Li4 1(a) 


0.1283 


0.0819 


-0.592 


Al4Li4 2(a) 


0.1276 


0.024 


-0.506 


Al4Na4 1(b) 


0.1302 


0.0909 


-0.174 


Al4Na4 2(b) 


0.1103 


0.0607 


-0.127 


AI4K4 1(c) 


0.0656 


0.0663 


-0.634 


AI4K4 2(c) 


0.0649 


0.0867 


-0.618 


C4H4 


0.245 


0.2410 


-0.030 


CqRq 


0.000 


0.2588 


-0.009 
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Table 2: The ground state dipole moment, /Ug, linear polarizability, a, 1st 
hyperpolarizability, {3 and the 2nd hyperpolarizability, 7, (tumbling aver- 
age) for the clusters and for trans- polyacetylene chain from ZINDO-MRDCI 
calculations. The units are in au. 'n' is the number of -CH — CH- units. 



Molecule 




a 




7 


Al4Li4 1(a) 


0.000 


4.9 X 10=^ 


542.5 


1.91 X 10^ 


AliLii 2(a) 


0.000 


5.5 X 10'" 


211.9 


5.33 X 10''^ 


Al4Na4 1(b) 


0.076 


5.9 X 10-^ 


8465.2 


1.09 X 10^ 


Al4Na4 2(b) 


8.6 X 10-^ 


8.7 X lO-'^ 


1098.5 


2.00 X 10^ 


AI4K4 1(c) 


0.004 


5.4 X lO-'^ 


79.3 


2.60 X 10^ 


AI4K4 2(c) 


5.720 


4.7 X 10^ 


1.2 X 10^ 


1.90 X 10^ 


C4H4 


0.000 


2.9 X 10^ 


0.000 


4.76 X 10^ 




0.000 


5.4 X 10^ 


0.000 


8.44 X 10=^ 


(CH=CH)„, n=l 


0.000 


136.3 


0.000 


2.78 X 10^ 


(CH=CH)„, n=2 


0.000 


421.0 


0.000 


4.15 X 10^ 


(CH=CH)„, n=3 


0.000 


852.4 


0.000 


6.17 X 10^ 


(CH=CH),, n=4 


0.000 


1455.2 


0.000 


2.82 X 10^ 


(CH=CH)„, n=5 


0.000 


2203.2 


0.000 


8.41 X 10^ 


(CH=CH)„, n=6 


0.000 


3074.9 


0.000 


2.07 X 10^ 
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Figure 1: Equilibrium ground state geometries for Al4Li4, Al4Na4 and AI4K4. 
The footnote of each structure contains the ground state energies in au. 
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Figure 2: Equilibrium ground state geometries for the other set of Al4Li4, 
Al4Na4 and AI4K4, very close in energy to Figure 1. The footnote of each 
structure contains the ground state energies in au. 
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